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Abstract 

Among the light nuclear clusters the a-particle is by far the strongest bound 
system and therefore expected to play a significant role in the dynamics of nu- 
clei and the phases of nuclear matter. To systematically study the properties 
of the a-particle we have derived an effective four-body equation of the Alt- 
Grassberger-Sandhas (AGS) type that includes the dominant medium effects, 
i.e. self energy corrections and Pauli-blocking in a consistent way. The equa- 
tion is solved utilizing the energy dependent pole expansion for the subsystem 
amplitudes. We find that the Mott transition of an a-particle at rest differs 
from that expected from perturbation theory and occurs at approximately 
1/10 of nuclear matter densities. 
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I. INTRODUCTION 

The modification of few-body properties such as the binding energy and the wave func- 
tion of a bound state due to a medium of finite temperature and density is an important 
subject of many-particle theory. As an example we consider symmetric nuclear matter con- 
sisting of nucleons (equal number of protons and neutrons) at density p and temperature 
T. The modification of singie-nucleon properties can be obtained from a Dyson equation 
in terms of a self-energy. In a specific approximation the quasiparticle picture can be de- 
rived. A more rigorous description leads to the nucleon spectral function. Similarly we 
can consider the two-nucleon system where the medium modification are obtained from a 
Bethe-Goldstone equation In addition to the self-energy shift, also Pauli-blocking has 
to be taken into account that is of the same order of magnitude. It has been shown ]2]|| 
that, as a consequence, deuterons in nuclear matter become unbound if the density exceeds 
a certain value, the Mott density. 

Of course, the same mechanisms are also responsible for the modification of higher clus- 
ters embedded in nuclear matter. However, the solution of the few-body in-medium equa- 



1 



tions where the effects of the medium are accounted for by a density and temperature 
dependent contribution to the Hamiltonian has only been done within perturbation theory, 
see 0. Therefore the results achieved for the energy shifts and the Mott densities are only 
approximations. 

Recently, rigorous methods have been used to find solutions for the three-body prob- 
lem in nuclear matter [0-0. The Faddeev equations are extended to include the effects of 
the medium, and the corresponding Alt-Grassberger-Sandhas (AGS) equations have been 
solved [jlCfl . Different properties of the three-nucleon system in the medium such as the 
modification of the binding energy of the three-nucleon bound state || and the medium 
modification of the nucleon-deuteron break-up cross section |5|||] have been calculated. 

In the present letter we give first results of the solution of the in-medium four-particle 
equation describing the modification of the binding energy of the a-particle in symmetric 
nuclear matter. An AGS-type equation has been solved and the results will be compared 
with those of perturbation theory. 

Note that the four-particle correlations in low-density nuclear matter are very important 
because of the large binding energy of the a-particle. They have to be accounted for not only 
in equilibrium when considering the nuclear matter equation of state or the contributions of 
correlations to the single-nucleon spectral function, but also in nonequilibrium such as the 
light cluster formation in heavy ion collisions. 



II. IN-MEDIUM FEW-BODY EQUATIONS 

The few-nucleon problem in nuclear matter can be treated using Green function ap- 
proaches. Within the cluster- mean field expansion [|j], a self-consistent system of equations 
can be derived describing a n-nucleon cluster moving in a mean field produced by the equi- 
librium mixture of clusters with arbitrary nucleon number m. A Dyson equation approach 
to describe clusters at finite temperatures and densities has been given in [|TTJ. However, the 



self-consistent determination of the composition of the medium is a very challenging task 
that is not solved until now. We will perform the approximation where the correlations in 
the medium are neglected so that the embedding nuclear matter is described by the equi- 
librium distribution of quasiparticles (see also || for the two-particle problem, or |5|-|8j for 
the three-particle problem). 

The extension of this formalism to describe n-nucleon correlations in nuclear matter 
will be given elsewhere. Here, we will give some of the basic relations, which are direct 
generalizations of the three-particle case. 

Let the Hamiltonian of the system be given by 

H = J2 §^ a ± ai + \ Yl v ^ u > 1/2 ') <A<&«v*v (!) 
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where a± etc. denotes the Heisenberg operator of the particle that includes quantum numbers 
such as spin s\ and momentum k\. The free resolvent Go for an n-particle cluster is given 
in Matsubara-Fourier representation by 

G (z) = (z - Ho)' 1 N = Bo{z)N, (2) 
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where Go, Hq, and N is a compact notation for matrices in the space of n particles with 
respect to the particle indices given below. Here z denotes the Matsubara frequencies z\ = 
7rA/ (—i(3) + \i with A = 0, ±2, ±4, ... for bosons and A = ±1, ±3, ... in the case of fermions. 
To simplify the notation we have further dropped the index n on the matrices, however use 
it if explicitly needed. The effective in-medium Hamiltonian Hq for noninteracting quasi- 
particles is given by 

n , 2 n 

^ = £^ + s * = 5> (3) 

i=l i=l 

where the energy shift and the Fermi function fi are 

Si =£V 2 (12,f2)/ 2 , (4) 

2 

^ = f( £ l) = e{£l -jk B T +1 - ( 5 ) 

The notation 12 means antisymmetrisation. The factor N in Eq. (fj) resembles the Pauli 
blocking or normalization of the Green functions. This factor is different for the different 
clusters considered depending on the number of particles n. It is given by 

N = hh...f n ±hh...f n (6) 

where / = 1 — /. The upper sign is for an odd number of fermions (Fermi type) and the 
lower for an even number of fermions (Bose type). Note that NR = R N. 

The full resolvent after Matsubara- Fourier transformation may be written in the following 

way 

G(z) = (z-H -V)- 1 N = R(z)N J (7) 
where the potential V is a sum of two-body interactions between pairs a, i.e. 

^ = E y «=E J w, (8) 

a a 

and V 2 a is the two-body potential given in Eq. (|l]). The sum runs over all unique pairs in 
the cluster. Note, that as a consequence of Eq. (||) ^ V, also R(z)N ^ NR(z) that later 
on leads to right and left eigen-vectors. 

To be more specific: If the interaction is between particle 1 and 2 (in the pair a = (12) 
of a cluster of n particles) the effective potential of Eq. (§) reads 

(12|Af 2 V 2 (12) |l'2'> = (/ 1 / 2 -/ 1 / 2 )V 2 (12,l , 2'). (9) 

A useful notion is the channel resolvent G a (z) for an n particle cluster, where only the pair 
interaction in channel a is considered. This may be written as 

G a (z) = (z-H - V^N (10) 
= ( z - H - N^V^N = R a (z)N. 
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Using Rq 1 (z), R~ 1 (z), and R~ 1 (z) it is possible to formally derive the resolvent equations 
in the standard way. To keep the formal equivalence to the isolated case, the n-particle 
channel t-matrix T a is defined by 

R a (z) = Rq{z) + R {z)T a {z)Ro{z) . (11) 

With the use of T a (z) = N 2 T 2 (z) Eq. (|TT| ) leads to the well-known Bethe-Goldstone equa- 
tion i 

T?(z) = V 2 a + V 2 a R (z)N«T°(z) = V 2 a + T 2 a (z)R (z)N«V 2 a . (12) 

We remark that similiar equations have been written down by various authors previ- 
ously[Q,0. 

Note that the above equations are also valid for the two-particle subsystem embedded 
in a larger cluster (three, four, or more particles). As for the isolated equations the effects 
of the other particles appear only in the Matsubara frequencies z (energies) of the other 
particles in the cluster. No additional blocking factors N related to the larger cluster arise. 
Also note, that the changes due to the Pauli blocking are in the resolvent Go not in the 
potential V2. However, it is possible to rewrite this equation and introduce an effective 
potential as seen in Eq. ([12] ) and use unchanged resolvents instead. Making use of the more 
intuitive picture of a blocking in the propagation of the particles (related to the resolvents) 
we find by Eq. ( |T2|) the correct expression for the t-matrix that enters into the Boltzmann 
collision integral (see, for example, [1|]). 



The derivation of the three-body equation is straight forward and has been given else- 
where 0-0. The AGS operator Up a (z) |10| for the three particle system is defined by 



R(z) = 5 Pa R a {z) + Rp{z)U p<x {z)R a {z). (13) 
Inserting Eqs. (|TI|) and ( |T2"D in the above identity we result with the AGS-type equation 

U Pa {z) = S^Roiz)' 1 + Y,^N2T2(z)R (z)U^(z), (14) 

7 

that includes now medium effects as Pauli blocking and self energy shifts. We used the 
notation 5 a p = 1 — 8 a p. This equation solves the three-body transition operator for a three- 
particle cluster as well as for a three-particle cluster embedded in a more-particle cluster, 
i.e. the effect of the other particles in the cluster is again only in the Matsubara frequency 
(energy) z. The definition of the transition operator given by Eq. ( ]T3"| ) was chosen so that no 
additional factor N appears in the final equation. This guarantees that the cluster equations 
are valid also if they are part of a larger cluster. Thus, the two-body subsystem t-matrix 
entering in Eq. ( |14"|) is the same as the one given in Eq. (0). Therefore, it is possible to use 
all results of the few-body 'algebra', in particular those based on cluster decomposition. 

The in-medium bound state equation for an n-particle cluster follows from the homoge- 
neous Lippmann-Schwinger equation and is given by 

|<M = Ro(E B )V\i> B ) = R^Y.N^^b) (15) 
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where the sum is over all unique pairs in the cluster. As shown in Ref. |8j for the three-body 
bound state it is convenient to introduce form factors 



\F P ) = Y,~h^V^ Ba ) (16) 

7=1 

that leads to the homogeneous in-medium AGS-type equation 

3 

\F a ) = Y,&«P N i T iMBz)\Fp). (17) 

0=1 

We may generalize the AGS method given in Refs. [|T^-[I7| to the in-medium four-body 

case 
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|^> = Rv{B A )N^{B A )R,{B A ) '(B^), (3 = 1, . . . ,6 (18) 

7=1 

where 



\<l)p)=R Q (Bi)Niv£\il)B A ). (19) 

Introducing the 3+1 and 2 + 2 cluster decomposition of the four-body system, denoted 
by t, a, . . ., the sum on the right hand side of Eq. fll8|) may be rearranged by introducing 
four-body form factors 

r a 

with (3 C cr, <5^ Q = ^q,, if /3, a C r and <5^ a = otherwise. The homogeneous in-medium 
AGS-type equation for the four-body form factors is then written 

TO = ^^(B^RoiB^T^iB^RoiB^), [3 c <r, 7 C r. (21) 

T7 

The driving kernel consists of the in-medium two-body t-matrix defined by the Bethe- 
Goldstone equation and the in-medium AGS-type transition operator defined in Eq. (|H|). 
Note also that an additional Pauli blocking factor N% occurs. 

The equations for the three-body scattering and bound state problem have been solved 
numerically in Refs. ||,|||8]|| . An exploratory calculation to study a possible a-like conden- 
sate (quartetting) has been carried out using a variational ansatz for the (2+2) channel and 
by neglecting the (3+1) channel ||19|| . 

Because of the medium dependence of the equations the calculation time increases dras- 
tically. This is due to the fact that the positions of the deuteron pole as well as the three- 
nucleon pole vary with the intrinsic momentum and are not fixed at the usual binding energy 
because of the phase space occupation through other particles. Presently a sufficiently fast 
and accurate method to solve the three- and four-body equations relies on the separabil- 
ity of the subamplitutes that appear in the AGS equations. To solve the four-body bound 
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states we utilize the energy dependent pole expansion (EDPE) |L8| that needs to be adjusted 
to the in-medium case, because of different right and left expansion functions due to the 
nonsymmetric effective potential. 

To be more specific, we assume the following expansions for the amplitudes of the re- 
spective sub-systems embedded in the four-body equation. For the two-body sub-system we 
have 

T -y( z ) - E \^-yn(z))t^ n (z){T yn (z)\ - E \9jn)t jn (z)(g in \ = ^ N^\g in )t in (z) {g in \. (22) 

n n n 

The last equation of the right hand side is used in the present calculation and reflects a 
simple Yamaguchi ansatz for the form factors [20]. The Pauli blocking factor then appears 
explicitly. This has been used in Refs. [0-0,0 and a comparison to the Paris potential is 
given in Ref. ||. For the present purpose this approximation seems sufficient. For the 
three-body subamplitudes we use the EDPE expansion 

t,/MU 

with 

|f#M> = (gan\Ro(z)\~gfSm)t Pm (B 3 )\t^) . (24) 

The Sturmian functions corresponding to the fixed energy B3 are given by 

f;!f> = J2(9an\Ro(B3)\gprn)tp m (B 3 )\tZ) (25) 



Vt. 



Vt. 



A' I 



(3m 



r^> = ^(9an\MB 3 )\g 0m )t Pm (B 3 )\T^) . (26) 

(3m 



Introducing the form factors 

l F "> = E^(^)IW( 5 4)(^(5 4 )|i2Q(S 4 )|^) (27) 

/3m 

we obtain the following homogeneous system of integral equations 

K) -EEE^^rwivwif^^)) t%{B,) if;'). (28) 

rt vk ■yn 

Formally these equations resemble the structure of the isolated four-body equations. How- 
ever, the dominant features of the influence of the medium, i.e. the self-energy correction 
and the Pauli blocking, are systematically taken into account. 

Inclusion of spin-isospin degrees of freedom and symmetrization is a challenging task for 
the four-body problem and done as for the isolated case. To this end we have intoduced 
angle averaged Pauli factors as explained e.g. in Ref. and fit the self-energy by use of 
effective masses. 
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III. RESULTS AND CONCLUSION 



To solve the four-body equation numerically we use a Yamaguchi type rank one potential 
for the 3 Si and 1 Sq channels. The parameters are taken from an early work of Gibson and 



Lehman pl| . We renormalized the calculated binding energy of the a-particle so that it 
coincides for the isolated particle with the experimental one. Presently, instead of using a 
more elaborated approach to the isolated four-nucleon problem we merely study the change 
of the binding energy due to the density and temperature of the surrounding nuclear matter. 
From our recent results for the three-body systems || , we argue that the change is not very 
sensitive to the particular form of the potential. We shall, therefore, leave the study of 
model dependences for a future communication. 

We calculated the binding energy of an a-like cluster with zero center of mass momentum 
in symmetric nuclear matter at temperature T = 10 MeV as a function of the nucleon density. 
The results are shown as a solid line in Fig. 1. The Mott transition occurs at a single particle 
density of pMott = 0.0207 fm -3 . For comparison we have given a perturbative calculation 
shown as dashed line. This calculation is based on a simple Gaussian wave function for the 
a-particle with the width fitted to the electric rms radius. Also the binding energy has been 
renormalized to the experimental value. The Mott density in this case is at 0.0305 fm -3 
that strongly differs from the value gained from the solution of Eq. (^Hj) . 

The corresponding curves for the triton and the deuteron are shown as dotted and dashed 
dotted lines respectively. Note, that these binding energies are for clusters at rest in the 
medium. Although it is an interesting case, when the sub-clusters embedded in the larger 
cluster vanishes as a bound state the sub-clusters in this case have a dynamical binding 
energy that depends on the cm. momentum. Nevertheless, the question of Boromenian 
states and the Effimov effect needs further investigation. 

Unlike for the triton the a-particle still exists at densities where the Pauli blocking factor, 
see Eq. (H), becomes negative. The usual procedure of symmetrizing the effective potential 
by proper square root factors fails. Therefore when solving the four-body equation we have 
to keep track of the right and left eigenvectors in the subsystem. 

In conclusion, we derived and solved for the first time an effective in-medium four- 
particle equation of the AGS type. Applying it to symmetric nuclear matter, we found that 
the binding energy of the a-particle decreases with increasing density due to Pauli blocking 
and disappears at a critical value of the density (about 1/10 of the nuclear matter density for 
T = 10 MeV). The dependence of the results on temperature and center of mass momentum 
will be the subject of an extended work. 
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FIGURE CAPTIONS 



Fig. 1. Binding energy of an a-like cluster with zero center of mass momentum embedded 
in symmetric nuclear matter at a temperature of T = lOMeV as a function of nucleon 
density. Solid line: Yamuguchi potential, renormalized to experimental binding energy at 
zero density. Dashed line: perturbation approach. For comparison, the medium dependent 
binding energies of the deuteron (dashed-dotted) and triton (dotted) are also shown. 
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FIG. 1. Binding energy of an a-like cluster with zero cm. momentum embedded in symmetric 
nuclear matter at a temperature of T = lOMeV as a function of nucleon density. Solid line: 
Yamuguchi potential, renormalized to experimental binding energy at zero density. Dashed line: 
perturbation approach. For comparison, the medium dependent binding energies of the deuteron 
(dashed-dotted) and triton (dotted) are also shown. 
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